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Abstract 
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■^j- , The exchange of light pseudoscalars between fermions leads to a spin- 

(N, 

independent potential in order g 4 , where g is the Yukawa pseudoscalar-fermion 

OS 
ON 



coupling constant. This potential gives rise to detectable violations of both 

the weak equivalence principle (WEP) and the gravitational inverse-square 

Qh! law (ISL), even if g is quite small. We show that when previously derived 

WEP constraints are combined with those arising from ISL tests, a direct 



experimental limit on the Yukawa coupling of light pseudoscalars to neutrons 

can be inferred for the first time (g^/An < 1.6 x 10 7 ), along with a new (and 

significantly improved) limit on the coupling of light pseudoscalars to protons. 
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In a previous paper ||T[, it was shown that laboratory bounds on the Yukawa couplings 
of light pseudoscalars to protons and neutrons could be significantly improved by using the 
results from recent weak equivalence principle (WEP) experiments [[|. These experiments 
are sensitive to the spin-independent long-range forces that arise in order g A from two- 
pseudoscalar exchange |j],[^||, where g is defined by the coupling 

C(x) = ig4>(x)j 5 4>(x)4>(x). (1) 

Here 4>{x) is the field operator for a pseudoscalar of mass m, and ip(x) denotes either a 
proton (p), electron (e), or neutron (n) of mass M p , M e , or M n respectively. For each pair 
of interacting particles, C(x) leads to a potential V^(r) in order g 4 which in the m — > 
limit is given by 0|| 

V^(r) = 9a9b (2) 

ab K } 64tt 3 M M 6 r 3 ' y } 

where a and b may each denote p, e, or n. The object of the present paper is to demonstrate 
that already-existing data from tests of the gravitational inverse-square law (ISL) ||,|] pro- 
vide new stringent constraints on gi and g\. When combined with the constraints implied 
by Eq. (^) below and the data from the WEP test in Ref. 0, the ISL data lead to the 
first direct experimental bound on the pseudoscalar- neutron coupling constant g^, and to a 
significantly improved bound on gi (see Eq. ( pp below). 

Leaving aside for the moment the contribution from electrons, it was shown in Ref. |T[ 
that Va b leads to an acceleration difference Aa2_2' of macroscopic test objects 2 and 2' in 
the presence of a common source M\. If these have masses M 2 and My, and contain Z 2 
(N 2 ) protons (neutrons), and Z 2 > (Np) protons (neutrons) respectively, then 



Aa 2 _ 2 ' = s{r) 



m 



h , 



m<± 



\/ i /2-2' \ r 1 I 2-2' J 



(3) 



In Eq. ([JD, F(f) is the integral over the mass distribution of the source |1|, /Xj = Mj/mn, 
m H = mdH 1 ), M n ~ M p = M, and A(Z//z) 2 _2' = Z 2 /fi 2 — Z 2 > / ji 2 i, etc., 0. Since all the 
parameters appearing in Eq. (Rt) are known, except for the pseudoscalar couplings gi and g\ 



(to protons and neutrons respectively), an experimental determination of Ac?2-2' leads to a 
constraint on gi and g\. 

As noted in Ref. |1|, however, the right-hand side of Eq. (|3]) vanishes whenever gi and 
9n satisfy 

gj_ A(iV///) 2 _ 2 , 

in which case #? and g\ can be arbitrarily large and still be compatible with any experimental 
bound on Aa2-2' • Since the right-hand side of Eq. (Q) is close to 1 for most pairs of materials, 
including those used in Ref. || , Eq. @ can be satisfied even when gi and g\ are each quite 
large provided gi ~ g 2 n . This is shown graphically in Fig. 1 which plots the constraints 
in the gi-g\ plane that emerge when Eq. (|3]) is combined with the experimental limits of 
Gundlach, et al. 0. It is seen that the boundary of the allowed region is a hyperbola 
with an asymptote near gi = g 2 , along which no limits on gi or g\ can be inferred. To 
circumvent the problem caused by such "hyperbolic" constraints, one can combine results 
from experiments using different materials, which thus have slightly different asymptotes. 
Alternatively one can choose special materials (such as 2 = Li and 2' = Ru) for which Eq. (|J) 
can never hold |TJ], and which thus lead to ellipses in the gi-g\ plane. The combination of 
"elliptical" and "hyperbolic" constraints would then lead to separate bounds on g% and g\. 
As we now demonstrate, existing data from ISL tests also provide elliptical constraints on 
gi and g\ and, when combined with earlier WEP results, lead directly to the bounds quoted 
in Eq. (|T0|) below. 

Consider the ISL experiment of Spero, et al. |5||| in which a cylindrical Cu test mass 
is suspended at the end of a torsion fiber inside a larger hollow stainless steel cylinder. It 
can be shown that for infinitely long cylinders the Cu test mass will experience a force 
from the stainless steel cylinder only if the underlying interaction is not a pure 1/r 2 force. 
When the small (and calculable) end effects due to finite cylinders are taken into account, 
the experiment of Spero, et al. becomes a null test for the presence of new non-Newtonian 
inverse-power- law interactions, such as V^ b (r) in Eq. (§). We will use two convenient pa- 



rameterizations of power law potentials between two particles 1 and 2 f?||| : 

-*-(^)(?r- 

Here G^r is the Newtonian gravitational constant, r = 1 fm is an arbitrarily chosen length 
scale, Bi and B 2 are the baryon numbers for bodies 1 and 2 respectively, and A n and a n 
are dimensionless constants characterizing the strength of the interaction. When gravity is 
included, the total potential energy between these two point masses is given by 

G N M l M 2 



V tnt r) 



r 
G N M 1 M 2 



r \ n—l 

1 + a n 

r 



1+A, 



(6a) 
m P \ 2 [Bi\ (B 2 \ froY- 1 ] (6b) 



m H s \ViJ \fi2. 



where nip = JTlc/Gn is the Planck mass. The null results of the ISL test of Spero, et al. 
can then be used to set limits on a n or A n , after integrating the corresponding l/r n+1 force 
laws over the mass distributions of the Cu test mass and the stainless steel cylinder [BJ . The 
lcr limits implied by Spero, et al. for a n and A n are shown in Table I for several physically 
relevant values of n. The results in Table I, which were obtained by direct integration over 
the mass distributions of the interacting Cu and stainless steel cylinders, are in excellent 



agreement with those obtained previously by Mostepanenko and Sokolov JK| who used 
a phenomenological parameterization of the non- Newtonian interaction to constrain A n . 
Although we are specifically concerned with the case n — 3, other values of n are also 
interesting: n = 2 potentials can arise from 2-scalar exchange, as well as 2-photon exchange 
llfl , and n = 5 characterizes the 2-body potential from neutrino-antineutrino exchange [12| 



and the 2-pseudoscalar exchange potential with derivative coupling (which is applicable to 
axions) [|]. Note, however, that n = 1 is uninteresting since such a potential would not lead 
to a deviation from the inverse- square law, but only to a modified value of Gn (which would 
be difficult to detect). Table I also presents the lcr limits derived from the experiment of 
Mitrofanov and Ponomareva (MP) |13J which is a test of the ISL over the range 3.8-6.5 mm. 



In this experiment a modified Cavendish apparatus is used to measure the force between 
a mass A suspended at one end of a torsion balance, and a second mass B whose distance 
from A is varied. The experiment then compares the experimental value for the force ratio 
F(r{)/F(r2), where r% = 3.773(40) mm and r-i = 6.473(40) mm, to the calculated ratio 
expected assuming Newtonian gravity. We see from Table I that for the case n = 3, which is 
our concern in this paper, the limits implied by Spero, et al. are more stringent than those 
of MP, although for n = 4, 5, 6 the reverse is true. 

To extract constraints on the pseudoscalar coupling constants gi and g\ from Spero, 
et al, we begin by considering the interaction of two macroscopic objects separated by 
a distance r that is large compared to their dimensions. From Eq. (0) the total 2-body 
interaction energy V± 2 in order g 4 is given by 

vS>{r) = -sip h y> zi + #*) & z * + **•) • <7) 

which should be compared with the parameterizations of Eq. (Q) for n — 3. For the actual 
geometry of the Spero experiment one must integrate over the mass distributions of the 
inner test mass and the outer cylinder |J, so that 1/r 3 in Eqs. (|7j) and (^) is replaced by the 
appropriate average (1/r 3 ). By combining Eqs. (|7j) and (|5p for n = 3 the constraint implied 
by Spero, et al. can be expressed in the form 

L 2 - + 9 2 n-) U~ + 9l-) = ^G N MV H a,rl (8) 

Using the la limits from Spero, et al. presented in Table I we then find 

(0.469^ + 0.540(£) (0.460^ + 0.549^) < 3.5 x 10" 12 . (9) 

In Eq. ([|) the expression in the first set of parentheses arises from the hollow cylinder, where 
the coefficients of g^ and g\ are the values of Zj \i and N//I, respectively, for stainless steel. 
Similarly, the expression in the second set of parentheses represents the contributions from 
the Cu test mass, while the right hand side of Eq. (|S|) is derived from the bound on a^r^ 
quoted in Table I. 



We see immediately from Eq. (§) that the constraint implied by the ISL experiment 
of Spero, et al. leads to an ellipse in the gi-g\ plane, as can be seen in Fig. 1. This 
is, of course, related to the fact that the left side of Eq. (|9]) cannot vanish unless both 
gi and g\ do. Figure 1 also exhibits the previously derived WEP constraint jl] from the 
experiment of Gundlach, et al. |fj], which gives rise to a hyperbola in the gi-g\ plane as we 
have noted previously. The significant new feature of Fig. 1 is that the combination of the 
hyperbolic WEP constraint and the elliptical ISL constraint lead to upper bounds on gi and 
g\ separately. We find from the figure the following la limits: 

g 2 p /4n < 1.6 x 10~ 7 , (10a) 

gl/An < 1.6 x 1(T 7 . (10b) 

The result for g 2 n in Eq. ([H]) represents the first direct laboratory constraint on the Yukawa 
coupling of pseudoscalars to neutrons. We note that the only previous laboratory limit 
on g 2 [[I] was based on an indirect model-dependent argument due to Daniels and Ni |TJ 



utilizing the spin-dependent results of Ritter, et al. JHJ. For gi there is an earlier result due 
to Ramsey fl6| , gi/An ^S 2.5 x 10~ 5 (lcr), which was obtained from a study of the molecular 



spectrum of H2. As we see from Eq. (^0|), the limit implied by combining the ISL and WEP 
results improves the Ramsey limit by more than two orders of magnitude. 

We can also extract from Fig. 1 constraints on gi and g 2 in special cases of interest. For 
a universal coupling to baryon number gi = g 2 , and we find (at the lcr level) 

glJA-n < 1.5 x 10- 7 . (11) 



The similarity of the results of Eqs. ( |T0"D and ( [LTD arises because the largest allowed values of 
gi and g 2 lie near the line g 2 = g 2 l1 as can be seen in Fig. 1. Two other results of interest are 
the limiting cases g 2 3> g\ and g\^> g 2 . We can see from Fig. 1 that these are determined 
by the WEP results of Gundlach, et al., and hence can be taken over from our previous 
analysis Jl|: 

gl/A-K < 9 x HT 8 , {gl » g 2 n ), (12a) 



gl/A-K < 7 x lO" 8 , (<£>$. (12b) 

Although we have focused thus far on the pseudoscalar couplings to protons and neutrons, 
it is straightforward to show that the contributions from electrons can be incorporated via 
the substitution g 2 — ► g 2 = g 2 + (M/M e )g 2 , where g e is the pseudoscalar-electron coupling 
constant. This follows from Eq. (Q) by noting that the contributions from electrons are 
enhanced by a factor M/M e relative to those from protons and neutrons. The limits on gij^K 
in Eqs. (|8|)-(|l2l) can then be taken over immediately for g 2 /4n, and these lead to constraints 
on g 2 , at least in principle. In practice, however, existing limits on g 2 /^ obtained from 
spin-dependent experiments Jl7j are more stringent, g 2 /4:ir < 10~ 16 . It follows that despite 



the enhancement arising from the factor M/M e , the contribution from the term in g 2 /4:7r 
proportional to g 2 is at most of order 1CT 13 . Thus, the bounds on g 2 /4:7i implied by Eqs. (§)- 
(|12|) are in fact bounds on g 2 /4:ir, and the prospects for constraining g 2 via V^ b seem quite 
remote at present. 

The limits on g 2 and g 2 in Eqs. fl9D-(|T2"D are the most restrictive direct laboratory con- 
straints currently available. Although astrophysical arguments based on stellar cooling cal- 
culations are more stringent [|H|, typically g 2 /4:7r < 10~ 21 , they are necessarily more model 
dependent. For derivative-coupled pseudoscalars such as axions, there is at present no viable 
alternative to astrophysical bounds, since those arising from existing laboratory experiments 
(corresponding to n = 5 in Table I) are too weak to be of use. However, by adapting the 
present formalism future laboratory experiments carried out over shorter distance scales may 
give rise to useful bounds on axions, as we will discuss in more detail elsewhere. 

In summary, we have shown that the limits implied by the ISL experiment of Spero, et al. 
0|, complement those previously derived from the WEP experiment of Gundlach, et al. [0, 
and together allow the pseudoscalar-neutron coupling constant g 2 to be directly determined 
for the first time. In addition, the combination of these two experiments leads to a new bound 
on the pseudoscalar-proton coupling constant g 2 , which improves on the earlier Ramsey limit 
fl6|j by more than two orders of magnitude. As was noted in Ref. |1[, the Gundlach results 



lead to hyperbolic constraints on gi and g\ which admit the possibility that each of these 
could be quite large, provided that gi ~ g\. Absolute bounds on these constants could be 
obtained, however, by using special materials such as Li and Ru. What we have shown here 
is that already existing data from the ISL experiment of Spero, et al. provide the needed 
elliptical constraints, in effect playing the same role that a WEP experiment utilizing Li 
and Ru would. The fact that the WEP and ISL experiments complement each other in this 
way raises the possibility of a new generation of WEP and ISL experiments whose results, 
when combined, would lead to even more stringent constraints on gi and g\. Although the 
laboratory constraints may not be as restrictive as those implied by astrophysical limits |18| . 



they are completely model- independent, and furthermore allow gi and g\ to be separately 
determined. 

The authors are deeply indebted to Riley Newman for providing the limits on a n rg _1 
quoted in Table I, and for numerous helpful discussions. This work was supported in part 
by the U.S. Department of Energy under Contract No. DE-AC 02-76ER01428. 
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TABLES 



Spcro, et al. 



A,-, 



Mitrofanov and Ponomareva 



An 



1.3 xlO -6 m 
1.3 xl0~ 8 m 2 
1.7 xlO -10 m 3 
2.3 xlO~ 12 m 4 
3.2 xlO~ 14 m 5 



7.7 xlO -30 

7.7 xlO~ 17 
9.9 xlO" 4 
1.4 xlO 10 

1.8 xlO 23 



6.8 xlO -5 m 
8.1 xlCT 8 m 2 

1.3 xlCT 10 m 3 
2.1 xlO" 13 m 4 

3.4 xlO" 16 m 5 



4.0 xicr 28 

4.7 xlO" 1(i 
7.5 xlO" 4 
1.2 xlO 9 
2.0 xlO 21 



TABLE I. la limits on a n r^ and A„, in Eq. (|J) from Spero, et al. [5] and Mitrofanov and 
Ponomareva [13]. 
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FIGURES 
FIG. 1. Laboratory constraints on gz and g\. The region in the g\—g\ plane above and to the 
right of each curve is excluded at the \a level by the indicated experiment. The gray shading indi- 
cates the region excluded by the overlap of all present laboratory experiments, and the remaining 
allowed region is shown in white. The data are from Gundlach, et al. [2], Ramsey [16], Ritter, et 
al. [15], and Spero, et al. [5]. The limit from Ritter, et al. on g\ is derived in Ref. [1]. 
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